Doubly excited ferromagnetic spin-chain as a pair of coupled kicked rotors 
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We show that the dynamics of a doubly-excited Heisenberg spin-chain, subject to short pulses 
from a parabolic magnetic field may be analyzed as a pair of quantum kicked rotors. By focusing on 
the two-magnon dynamics in the kicked XXZ model we investigate how the anisotropy parameter - 
which controls the strength of the magnon-magnon interaction - changes the nature of the coupling 
between the two "image" coupled Kicked Rotors. We investigate quantum state transfer possibilities 
and show that one may control whether the spin excitations are transmitted together, or separate 
from each other. 
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I. INTRODUCTION 

Over recent years, there has been sustained interest 
in coupled quantum systems. Numerous studies investi- 
gated the causes and effects of decoherence on a subsys- 
tem as it becomes entangled with its environment; others 
probed the generation of bipartite entanglement between 
a pair of quantum systems. It is vitally important to un- 
derstand these processes so they can be accounted for in 
protocols for quantum computation and communication. 

Studies of decoherence also shed light on the emer- 
gence of classical behavior from quantum dynamics P, Q ■ 
Quantum systems with a chaotic classical limit often fea- 
ture in such studies. For example, they can play the role 
of the environment: a ID system which displays chaos 
can replace a many-body heat bath (often modeled by 
an infinite collection of quantum Harmonic Oscillators) 
as a source of decoherence [3] . Other studies focused on 
entanglement generation: the rate of growth of the von 
Neumann entropy of the subsystem - ie the rate at which 
the subsystem becomes entangled with its environment 
- is directly related to measures of the chaos in the sub- 
system's classical limit 0, H, @1 ■ 

The chaos paradigm known as the Quantum Kicked 
Rotor (QKR) 3 plays a central role in these studies. The 
QKR corresponds to the dynamics of independent quan- 
tum particles evolving under the rather simple Hamil- 
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tonian Hi = ^ -\- K sniXi^^5{t — nT), where K rep- 
resents the kick-strength and T the kick-period. Cold 
atoms in pulsed standing waves of light were found to 
provide a very clean realization of the QKR: in 1995, 
the phenomenon of the quantum suppression of classi- 
cal chaotic diffusion was clearly demonstrated experimen- 
tally 3; 

later, the recovery of the classical diffusive be- 
havior in the presence of decoherence was also observed 
0. These works were followed by other studies by differ- 
ent experimental cold-atom groups worldwide 'l(J| prob- 
ing wide-ranging aspects of the QKR dynamics. In a pre- 
vious work [uJ, [iJl , we proposed that the singly-excited 
Heisenberg spin-chain in a pulsed parabolic field could 



provide an exact physical realization of the QKR: the 
dynamics of the spin- waves are given by a time-evolution 
operator of analogous form to that of the QKR. 

Coupled QKRs, have also been investigated in a num- 
ber of theoretical studies, though, unfortunately, no 
physical realization has yet been achieved. In this case, 
one considers two QKR Hamiltonians with an additional 
coupling potential V, i.e. H = Hi + H2 + F(a;i, 2:2, i). 
In 13], interactions which depended on the separation of 
two rotors with a non-local sinusoidal term were investi- 
gated; in [T3| the two particles were confined to within 
a short distance of each other. However, several stud- 
ies considered a sinusoidal coupling term dependent on 
a center-of-mass coordinate, [3 Ej [13 such as e.g.: 
V{xi,X2, t) = K12 cos (a:i -\- X2) En ^ " "-^)- 

In this work we show for the first time that the doubly- 
excited Heisenberg spin-chain system may-to a good 
approximation-be analyzed as pair of coupled kicked ro- 
tors. In fact, in this system. Nature even provides a cou- 
pling term of the centre-of-mass form K12 cos (xi -\- X2). 
The mapping is, however, far less straightforward than 
that found for the one-excitation system in [ll|, [l2| : the 
coupling here is mediated by bound-pair eigenstates (not 
found in the corresponding one-excitation chain), rather 
than spin waves, so acts only over a restricted part of the 
"image" phase-space. The wavenumbers of the bound 
states are complex, further complicating the mapping. 
Nevertheless, the analogy holds sufficiently well, so one 
can use it to shed insight on the dynamics. Further, it 
points to useful applications in state transfer, since we 
can use this understanding to control whether the two 
spin-flips propagate along the chain together, or sepa- 
rately. This adds to other applications that make use of 
the single-excitation correspondence [Tsj . 

In Section II we review briefly the one-particle dynam- 
ics of the Heisenberg XXZ spin chain and its mapping to 
the QKR. In section HI we consider the doubly-excited 
spin chain. We summarize essential features of the well- 
known XXZ eigenstates and their dependence on the 
anisotropy parameter A. We then introduce the anal- 
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ogy with the two-particle coupled QKR and explore the 
dynamics when the initial state consists of two neigh- 
boring spin-flips. We also highlight two cases where the 
kicked rotor correspondence simplifies: i) when A = 
the kicked spin-chain can be mapped to a pair of inde- 
pendent QKRs; and ii) when A ^ 1 the bound states 
effectively trap two excitations on neighboring sites and 
we show that in this limit, these bound states give rise 
to a further analogy with the QKR. We finish, in sec- 
tion IV with examples of how we can use these results to 
manipulate correlations in the spin- flip locations. 



II. THE HEISENBERG SPIN-CHAIN AND ITS 
ONE-PARTICLE IMAGE 

The well-known spin- 1/2 Heisenberg XXZ chain is gov- 
erned by the Hamiltonian: 

Hhc = -— ^ (2{a:^a~^;^ + anf^n+i) + Aaf crf+j^ af . 

Ti — 1 n 

(1) 

When investigating the dynamics of spin-chains such as 
this, a useful approach is to invoke quasi-particle models 
and interpret excited states as systems of indistinguish- 
able particles. In some cases, it is even possible to map 
the dynamics to a one-body "image" system which ap- 
proaches the classical limit as iV ^ oo [19|. 

iJhc conserves the number of spin-flips and a single ex- 
citation represents a spin- wave which distributes a single 
spin-flip throughout the chain. Higher excited-states cor- 
respond to multiple spin- waves which interact when they 
coincide through both an exclusion process (no two spin 
flips can simultaneously occupy the same site) as well as 
a mutual interaction induced by the a'^cr^j^-^ (Ising) term 
- the strength of which is determined by the anisotropy 
parameter A. Note that A = corresponds to the XXO 
chain and A = 1 is the isotropic Heisenberg chain. 

The eigenstates for the single spin-flip sector of ([T|), 
spanned by the basis states {\n) = ]]...]) : n — 
1, . . . ,N}, are translationally invariant magnon states 
with momenta k: 

i'^)-^E^'™i")- (2) 

Note that periodic boundary conditions are used, i.e. the 
configuration is a closed ring with cr°^_^ = a". The 
magnon momenta are determined by these conditions 
and take the values k = 2nI/N, I = 1,...,N. These 
states have energy: 

E-Eo = 2B + J{A-cosk), (3) 

where the ground state energy Eq = — JAiV/4, i.e. 
Hhc\ TT ■■• T) = Eq\ tt ••• T)- Adding an external 
parabolic kicking field to the Heisenberg Hamiltonian 



gives: 

B ^ 

H^Hu.^^Y.^n- n^fal ^ bit - jT), (4) 

n=l i=l 

where the kicking field has amplitude Bq with minimum 
at ng. Time evolving the time-periodic H for one period 
T yields a unitary map, 

m^{j + l)T)) = U{T)m=jT)) (5) 

where 

U{T) = e~'^^/'* = e*T^ T."=t{n-naf,yf,^-iTH,,Jh 

since the (5-kick nature of the time-dependent field per- 
mits us to split the operators. 

Using ^ and Q, it was shown in [lllj[l3|i that the ma- 
trix elements of U{T) in the single- flip basis {|n)} have a 
form very similar to the matrix used to evolve the quan- 
tum chaos paradigm, the QKR. These are given by: 

C/„„, ~ e^^("-"«''*"'-"J„,_„(JT) (7) 

for the spin-ring (an analogous form was given in [Tlj 
for open boundary conditions). Here, Jn{x) denotes a 
Bessel function of order n and, for convenience, we have 
set h — 1. 

We recall the form of the QKR Hamiltonian: 

i?QKfl = y-fccosi^(5(i-jT). (8) 

i=i 

with X e [0, 27r). In its classical limit, the dynamics is 
described by the famous Standard Map which is known 
to display a transition from integrability to chaos as the 
Stochasticity Parameter, K = kT, is increased. For 
K < 1 diffusion in momentum is blocked by invariant 
tori running through classical phase space. At large K 
phase space is almost completely chaotic and unbounded 
diffusion in momentum is typically seen. However, in cer- 
tain ranges of X w 2j'7r, j G Z, small transporting islands 
known as "Accelerator Modes" (AM) appear in classical 
phase space and give rise to anomalous diffusion. In the 
QKR, diffusion of momentum at large K is suppressed 
by quantum interference in a process known as Dynami- 
cal Localization [2l[ . We can express the QKR time 
propagator in a basis of plane waves |/) = exp(iZa;): 

{l\UQKRiT)\l') = e-'''^/^/-'j,^i (^^) . (9) 

Here r = hT is the rescaled effective Planck's constant. 

Comparing the above with ([7]) we see that the QKR 
and spin-chain propagators are of similar form, provided 
we identify JT and note that the kicking field 

Bq T plays the role of an effective Planck's constant. 
In effect, the spin-chain equivalent to the if cos a; term 
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in Hqkr arises from the dispersion relation of the spin- 
waves, ie re- writing ([l}: 



Hhc = ^ J(A - cosk)|k)(k| 



(10) 



So, to make the QKR — s- spin-wave mapping we also had 
to identify position (x) in the QKR, with momentum 
in the spin-chain (k); and momentum in the QKR with 
position (spin-site) in the chain. 

With the aid of this mapping we can identify the spin- 
wave equivalent of classical chaos phenomena such as Ac- 
celerator Modes (AM), transport on tori ^22], cantori or 
stable islands, and quantum chaos phenomena like Dy- 
namical Localization The classical transporting is- 
lands represented by the AM have evident potential ap- 
plicability in quantum state transfer so, below, we in- 
vestigate these in particular: they occur for K « 2j7r 
where j is an integer. In the classical image phase-space, 
they correspond eg to initial conditions located around 
(xqjPo) ~ (i'''/2,0) which "hop" in momentum each pe- 
riod such that a,t t — riT (n e Z): 



(a;„,Pn) ~ (7r/2, T'^nnj). 



(11) 



Quantum mechanically, if the effective Planck's constant 
is small enough, these islands can support Gaussian 
states that follow the classical trajectories - i.e. they 
"hop" in momentum every period 23']. Gaussian exci- 
tations were indeed seen in the one-flip spin-chain, 
provided the initial spin-flip occurs at a site near the min- 
imum of the magnetic kicking field and K — JTBq « 
2nj. The excitations were seen to hop approximately 
2tt/Bq sites each period, with little dispersion. 
We now consider the two-flip case. 



and also a relation between 9 and the quasi-momenta, 
the Bethe Ansatz Equation (BAE): 



1 ^ g'*(Ki+«;2) _ 2Ae*'*i 
I gi(Ki-i-K2) _ 2Ae''*2 



(15) 



Further restrictions are imposed by the periodic bound- 
ary conditions: 



Nki ^ 27rAi + 6*, Nn-i = SttAj 



(16) 



where the Bethe quantum numbers Ai < A2 are integers 
in the range Ai G {0, 1, . . . , — 1}. By solving the cou- 
pled system of equations in (|15p and (|16p . ki,2 and d can 
be obtained. Broadly speaking, these solutions fall into 
two groups depending on whether d has an imaginary 
component. The majority of the solutions of p5)) are 
real - these correspond states of two magnons which scat- 
ter off each other. For A = all the available solutions 
of are real (and equal to tt). When B is complex, 
the eigenstates correspond to bound states of two spin 
flips. The probability amplitudes of these states are at 
a maximum when the flips are on neighboring sites and 
they decay exponentially with the separation of the flips. 
While for any given A 7^ the widths of these states vary 
with the total momentum ki -I-K2, they become narrower 
as A increases. Crucially, for long chains [N 00) the 
energy of these states can be written [25| : 

£; - £^0 = 45 + JA - ^(1 + cos(ki + K2)). (17) 
when A > 0. 



B. Analogy with a pair of coupled kicked rotors 



III. TWO SPIN EXCITATIONS 

A. Bound-pair states and spin-waves 

Eigenstates in the double excitation sector are ex- 
pressed, via the Bethe ansatz, as pairs of spin waves [2^ : 

K2) = ^(ki, K2) ^ a("-l,f^2) |f^l,f^2)• (12) 

0<Jii<n2<-'V 

where A{k\^K2) is a normalization constant. |rti,rt2) de- 
notes a state with a spin- flip at sites n\ and rii. Bethe's 
ansatz for the amplitude is 

a(ni,n2) = e^(''i"i+«2"2+«/2) e*(«i"2+K2ni-e/2)^ (--^3^ 

The scattering phase 0(ki,K2) accounts for the interac- 
tion between the pair of spin-waves. On applying H^c in 
([T|) to these states and solving the eigenvalue equations, 
one obtains the dispersion relation 

£■ - i;o = 45 + J(2A - cosKi - COSK2) (14) 



The departure point for our analysis of the spin dy- 
namics as a system of coupled QKRs is the two-excitation 
spin-Hamiltonian, equivalent of (|10p : 

Hhc = -Ps I]^^ 2 J (A - cos(ki) - cos(k2)) Iki, K2){«:i, K2I 
-hi(£;o + 4B-h AJ) 

(18) 

for A > 0. Here is a projector on to the scattering- 
state component of Hilbert space, and Pf, onto the bound 
states. 

Comparing the above with the typical coupled QKR 
potential V{x\,x-2) — K\ cos x\+Ki cos a;2 -I- A'12 cos(xi -I- 
x-i) and identifying Xi and Ari2 ^ might suggest 
that the scattering states be interpreted as giving rise to 
a pair of kicked rotors; and that a coupling between these 
rotors arises due to the bound states. However, we note 
the important difference that the n,i for the scattering 
and bound states correspond to complementary portions 
of the "image" phase space. For the bound states, 
is complex, but (ki -|- K2) is real. In addition, we show 
below (in ([H])) that in fact, for large A, Ari2, i.e. the 
effective coupling is twice as large as suggested by (IT5|) . 
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Figure 1: Showing the production of "bound state" accelerator modes (AM2) which move slowly, and fast "scattering state" 
accelerator modes (AM). When A = only the AM are present, in this case the dynamics maps to two independent QKRs. 
With increasing A the AM2 become dominant. We plot the on-site magnetization for for two initially neighboring spin-flips 
\ip{0)) = |400,401). We have chosen here, the parameters Ks — 13, Bq = 0.1, no = 400 and T — 1 for chains of 800 spins, 
which are known to produce Gaussian excitations in a singly-excited chain. 



The parabolic kick will couple the eigenstates to 
each other (including coupling bound-pair and scattering 
states). As A increases, the overlap between the bound 
and scattering state energies decreases and the two bands 
separate for A > 2. This will suppress the coupling and 
imply that for large A, if the initial state has negligible 
overlap with the bound subspace, the dynamics will be 
essentially uncoupled. 



C. Evolution of |n, n + 1) initial states 

In this section, we explore the dynamics of an initial 
state prepared with two spin-flips localized on neighbor- 
ing sites near the center of the chain \tp{0)) = \no,no + l). 
Parameters corresponding to accelerator modes are used: 
JT = 130 and Bq = 1/10, so if - Att. Fig. [T] shows the 
resulting on-site magnetization (P^) — (i(l — cr^f )) of 
IV'(O)) after successive applications of the map ([6]) for 
A = 0, 1 and 2. 

When A = 0, a pair of hopping wavepackets is pro- 
duced. Each travels 27r/i?Q « 130 sites each period. This 
is consistent with the single-particle accelerator modes 
(see previous section). In contrast, when A = 1 or 2, 
there are two sets of hopping wavepackets. One pair 
(AM) evolve like those in the A = chain, while the 
other pair (AM2) hop approximately tt/ (ABq) sites each 
period. For A = 1 the AM2 wavepackets decay rapidly 
and by the third period (t = 3T) they are almost indis- 
tinguishable from the chaotic central portion. 

To get a more complete picture of the dynamics we 
plot, in Fig. [21 the two-site correlation function (Pn-^Pn^) 
for \ip{T)), allowing us to follow the relative positions of 
the spin-flips. We find that the AM2 wavepackets con- 
tain flips that travel together, this suggests they are sup- 
ported by the bound states. The AM wavepackets on 
the other hand appear in an anti-correlated portion of 
the wavefunction. 



We now consider these two different kinds of behaviour 
in more detail. 



D. A = and 'independent' QKRs 

For A = only the exclusion interaction is present 
between flips. The effects of this interaction are subtle 
and sensitive to the initial conditions. For certain cases, 
where the flips are initially separated by an odd number 
of sites, it has been shown to change the character of 
entanglement when the two excitations collide [1^. A 
separate study on the transfer of entangled states in a 
doubly excited XXO chain was carried out in . 

Here, we are interested what influence the exclusion in- 
teraction has on the QKR-like behaviour of excitations in 
the kicked spin-chain. The A = model can be mapped 
to a system of spinless fermions via the Jordan- Wigner 
transformation (see appendix for details). The number 
of fermions matches the number of spin-flips. The exclu- 
sion interaction is accounted for by the anti-commuting 
property of the fermionic operators. Consequently, the 
fermions are "free" (non-interacting) and can be evolved 
separately under the single-particle dynamics. The re- 
sult of this is that the kicked chain maps to a system of 
non-coupled QKRs. However, in the spin representation, 
the exclusion interaction is still relevant (the spin-flips 
do not evolve separately). To see this we make use of 
the Floquet operator in the two spin-flip basis (see ap- 
pendix) : 

(rii,n2|C/''^°(r)|mi,m2> = 

g-»-|^[(ni-"o)^ + {"2-io)^] j"l+"2-mi-"i2 

[Jni—mi {f^)Jn2~m2 iP) ~ Jni—m2 {f^)'Jn2~mi {P)] 7 (1^) 

where /3 — JT and N ^ oo. 

The effects of the exclusion interaction are not ac- 
tually seen in Fig. 1. For example, the on-site mag- 
netization after one period (i.e. for \^{t = T))) is 
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Figure 2: (color online) Spin-spin correlations corresponding to Fig. [T]at t = T. The two-site correlation function {P„-^P„^), 
equal to the probability of finding one flip on site ni and the other on n2, is shown. At A = the spins are anti-correlated in 
contrast to Fig. [T] which suggests the dynamics of uncoupled particles. A = 1 and 2 have an anti-correlated component (flips 
separate) as well as an additional component where the flips travel together. 



which is (P^) 



|7-/A=0 12 , 

I ^ ni ,n,no ,no-t-l I 



^^n9 <n 



rjA=0 12 
n,n2.no,no + l I ' 

This is the 



same as for two independent spin-flips initialized at sites 
no and no -I- 1. Using the free-fermion correspondence, it 
is straightforward to show that for all later times (P^) is 
exactly equivalent to the sum of expectations for a pair 
of uncoupled QKRs. 

The coupling induced by the exclusion interaction 
is, however, evident in Fig. 2, which plots the two- 
site correlations after the first period: (Pi^Pij) = 
\^nun2 no no+il"^ ■ cffcct, for this particular initial 

state, is to prevent the spin- flips from travelling together. 
The two site correlation is highest when the flips travel 
JT = 130 sites in opposite directions. If the flips were 
non interacting (i.e. allowed to co-exist on the same site) 
then it would be equally likely the flips would travel to- 
gether or apart. Different correlations are seen when the 
initial separation of the flips is changed. 

So when A = 0, where the Heisenberg chain eigen- 
states consist entirely of scattering states, the behavior 
of two spin-flips is like that of two Kicked Rotors except 
the flips build up correlations in their relative positions. 



fective Hamiltonian when A ^ 1. In this approximation 
the bound state amplitudes are: 



a(ni, n2) = (S„i, 



-le 



i{Ki+K,2)ni 



(20) 



and their dispersion relation remains unchanged from 
P?)) . Clearly, in center of mass coordinates, the bound 
states have the same form as a single magnon solu- 
tion. Naturally, this similarity extends to the dynam- 
ics of states on the nearest neighbor (NN) subspace, 
{|n,n-|- 1)}: Two initially neighboring spin-flips evolve 
together in approximately the same way as a lone flip in 
the single excitation basis {|n)}. We anticipate that for 
A > 1: 



{ni,n2\Uhcit)\m, to -f 1) w i'' 



Jt_ 
2A 



■Jni .112 — 1 

(21) 

where the propagation of the neighboring flips is slower 
than for a single flip - it is scaled by J/(2A) rather J. 

The influence of the kicking field on the NN subspace 
can be incorporated into (pT|) to give: 



E. Bound State QKRs for large A 

When A 7^ 0, the spin-flips can form bound pairs, and 
when A = 1 or 2, a pair of neighbouring spin- flips will 
overlap with both bound and scattering eigenstates. The 
additional features seen in the probability distributions 
when A = 1 and 2 are remnants of QKR-like behavior 
of bound states that appears in the limit of large A. In 
this limit, the bound states confine the flips to neighbor- 
ing sites. Santos and Dykman [2^ use a perturbation 
expansion in spin coupling strength J to produce an ef- 



(ni,n2|i7ff(T)|TO,m + l) « 5„,,„,_ie'^«("i"("»-2)) 

■ni-rn j ( JT \ 

(22) 

up to a global phase. Again, we see an analogy to a 
QKR, with stochasticity parameter Kf, = JTBq/A and 
effective Planck's constant Tf, = 2Bq. 

We expect the accuracy of this approximation to fall 
with decreasing A as the bound states become broader 
and are coupled more strongly to the scattering states 
by the kicking field. However, we show in Fig. [3] that 
even for A = 2, QKR-like behavior is still seen on the 
NN subspace for short times. 
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Figure 3: Showing the decay, over time, of the correspondence 
between the dynamics of nearest neighbor flips and a QKR. 
F is the fidelity of the time evolution of two spin-flips intially 
on neighboring sites \Tp{t = nT)) = [U(T)]"\100, 101) to the 
matrix elements (|2H) and (|22|l for various Bq and JT/A = 5, 
A = 2. 



Scattering State QKRs and Center of Mass 
Diffusion 



We now consider parameter ranges for which a sin- 
gle particle displays Dynamical Localization. Taking 
K = JTBq = 5.0 and Bq = 1 a lone flip ini- 
tially spreads diffusively but at long times this spread- 
ing saturates and the flip becomes exponentially localized 
(P,|) ~ exp{— 2|n — riinitl/-^} with a locahzation length 
L — [JTY /A. The diffusion time is usually increased for 
coupled kicked rotors, e.g. in a related study [2§| a pair of 
rotors coupled locally in momentum U5pi , p2 were shown 
to localize with a much greater L. 

In Fig. |3]we follow the center of mass spreading of two 
spin-flips using the second moment {{ni+n2 — 2noY Bq) . 
The flips are initialized 10 sites apart so in the limit of 
large A this state should overlap only with scattering 
states. The spin-distribution localizes for A = as ex- 
pected for an uncoupled QKR; however, for larger A, 
the diffusion is not halted, but slows down appreciably 
after the "break-time" at A = 0. This slower diffusion 
saturates and reaches a constant rate for A > 1. 

For large A, due to the large energy gap, the kicking 
field will not significantly couple the bound and scat- 
tering states so the quantum state is supported only by 
the scattering states for all time. The behavior of the 
diffusion however, does not reduce to that of uncoupled 
kicked rotors (as might be suggested by the dispersion 
relation in UHl)). This is because the scattering states 
do not exactly reduce to a pair of magnons-even in the 
large N limit where the corrections to the momenta k.i/2 
(see (HH)) vanish (i.e. 9/N 0). They are distorted by 
the Ising term and correspondingly the presence of the 
bound states. This can be seen by rearranging the Bethe 
equations in section A. For A 7^ the scattering state 
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Figure 4: Influence of the erf coupling on the growth of 
the 'center of mass' second moments for two flips initialized 
10 sites apart and parameters K = 5 and Bq = 1 on a chain 
of 200 spins. 



amplitudes can be written 

a{ni,n2) oc e'~^"^'^"^^ [sin {nr {ni — n2 + 1)) 
cos ^ 



■ sin {Kr (ni — n2)) 



(23) 



where = 27r (Ai -I- A2) /A^ and = 

(tt (Ai - A2) + 61) /iV. Notably, when A > 1, the 
scattering states typically have no overlap with the 
\n,n + 1) subspace (except, eg, for the Ai = A2 = 6' = 
state) as this is occupied by the bound states. 



IV. POSSIBILITIES FOR CONTROLLING THE 
EVOLUTION OF SPIN-FLIPS 

These results may be of interest in the context of 
quantum information and state transfer as they suggest 
possibilities for manipulating the evolution of spin-flips 
(and spin-correlations) in a Heisenberg chan. Clearly, 
the evolution of a two-particle state on the non-kicked 
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XXZ chain depends on A and the shape of the initial 
wavepacket. These two factors also determine the pro- 
portion of the wavefunction that is supported by the scat- 
tering/bound states. In the kicked chain, for large A, the 
scattering and bound states correspond to different QKR 
images: The bound (6) and scattering state (s) QKR im- 
age parameters are related via Kb = Kg/ A and Tf, = 2ts. 
So, by picking suitable values of JT and Bq we can se- 
lect which dynamical regimes the bound and scattering 
components correspond to. 

For example, one could halt the propagation of cither 
the bound or scattering state portion of the wavefunction 
and allow the rest to travel. A possible way to do this 
is to make use of resonances in the QKR. These occur 
for T — Airr where r is rational. For r ~ I (primary res- 
onance) ballistic spreading occurs in momentum for the 
QKR (position for the spin chain) and when r = 1/2 (an- 
tiresonance) diffusion in momentum can be suppressed. 
These two conditions could be achieved for the bound 
and scattering states respectively by setting ti, — An. 
This would lead to ballistic diffusion for initially neigh- 
boring flips and could prevent flips that are initially well 
separated from spreading. 



Making use of the transformation and setting A = 0, 
the kicked spin-chain Hamiltonian ^ becomes: 



H = 



-4[e, 



JV-1 
3 = 1 



1+C, 



+ ^E"i(j-jo)^(i 



2c]ci) 5 (t/T) . 



(A2) 

The transformed Hamiltonian has boundary terms that 
depend on whether the number of fcrmions, r, is odd or 
even; these arise from the periodic boundary conditions 



'N+l 



and crffj^i = erf 



We now calculate the result of time evolving over 
one period. Using the Heisenberg picture, we define 
c](T) = W{T)cj{0)U{T), where U(T) is the Floquet op- 
erator of eq. [6] The absence of any mutual interaction 
(i.e. terms of the form f (jij2, Ji, i2)cjiCj2C^/ c^, ) in (|A2[) 

implies the fermions are "free" and therefore the fermion 
operators can be time evolved using single particle states: 
W{T)c{iIj{0))U{T) = c\ip{-t)) where c'^{tp[-t)) creates 
the single particle state ip{—t) = {T)il;{0). This cor- 
responds to the single-flip basis, so Cj{T)\(f)) = U'^{T)\j). 
Using the matrix elements t/„n' in ^ gives, up to a 
global phase. 



CONCLUSIONS 



c]{T) 



-jK- , .- e'^«(^'-^°)'/2^i^'-Vy_,(-jr)4(0). (A3) 



We have investigated the dynamics of a pair of spin- 
flips on a periodically kicked Heisenberg chain, focusing 
on the roles of the scattering and bound eigenstates of 
the underlying time independent model. Analogies to 
coupled and independent rotor systems have been iden- 
tified and analysed. 

T. B. acknowledges support from the EPSRC. This 
work was partly supported by Grant-in- Aid by MEXT, 
Japan. 



Appendix A: TIME EVOLUTION FOR THE 
KICKED XXO MODEL 

Here, we show how the kicked XXO chain (obtained 
from (HI) by setting A = 0) maps to a system of in- 
dependent QKRs. This is done by first applying the 
Jordan- Wigner transformation [so'l , a non-local mapping 
of spin-flips on the chain to free fermions on a lattice. 
This transformation defines fermion creation and annihi- 
lation operators, and Cj respectively, in terms of spin 

operators so that af = (1 — 2clci) and 



n 



(l-2c]c,0 



(Al) 



The product of (1 — 2cjcj) terms accounts for the differ- 
ence between intcr-particle exchange statistics - negative 
for fermions and positive for spin-flips. The creation and 
annihilation operators obey the standard commutation 
relations: {ci,cj} = 5ij and {ci,Cj} = and are defined 
with respect to a vacuum state |</)) such that Cj\(p) — 0. 



The equivalence between this propagator and the time- 
evolution for the QKR is clear and the parameters cor- 
respond as before: JTBq — > K and Bq r. There- 
fore, a kicked fermion evolves in position in the same 
way as a QKR evolves in momentum i.e. j — > It. This 
multiplc-fcrmion correspondence is a direct extension of 
the single-flip analysis jTllj • 

Calculating two-particle correlation functions is now 
straighforward. For example, we find the matrix ele- 
ment of the Floquet operator in the two spin-flip basis 
|ni,?i2) = -4^41 1 



(ni,n2|t/(r)|mi,m2) 



Substituting (jA3p into this 



\Cr,^Cr,^U{T)cl,^c'^^^\(t)) (A4) 

0|c„.,c„,4(T)ct jr)|0>j 



(ni,n2|(7(T)|mi,m2> = ^ ^-^^[i^^i-iof H^^-oof] ^ 



«l+*2— "-1— "-2 



Jil—ni ( JT^ ,Jj^2—n2 ( JT^ X 



(A5) 



From Wick'stheorem ((/)|c„iC„2cJ^cJj0) = 5,ni^i^5m2,i2~ 
Smi,i2^m2,ii ^^'^ therefore, 

(ni,n2|l7(r)|mi,m2) = e-"-^[("i-^"''+("^-^"''l x 

{-JT)Jm2-,n{~JT)]. 
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Substituting Jj{—j3) = J-j{j3) into this yields ([T9)) . As a result, the corresponding QKR images will be cou- 

Finally, we note that when A ^ the Ising term, pled. 
Hzz = — JA Yln'^n'^n+i/'^^ iutroduces a mutual inter- 
action between the fermions: 

N 

ffzz = JA ^ c]c, - 1/4 - c]c]^,c,+^c, . (A7) 
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